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Abstract 

The nonlinear sigma-model and its generalization on 
iV-component spins, the O(N) model, are considered 
to describe thermodynamics of a quantum quasi-two- 
dimensional (quasi-2D) Heisenberg antiferromagnet. A 
comparison with standard spin-wave approaches is per- 
formed. The sublattice magnetization, Neel temperature 
and spin correlation function are calculated to first order 
of the 1 /TV-expansion. A description of crossover from a 
2D-like to 3D regime of sublattice magnetization temper- 
ature dependence is obtained. The values of the critical 
exponents derived are j3 = 0.36, rj = 0.09. An account 
of the corrections to the standard logarithmic term of 
the spin-wave theory modifies considerably the value of 
the Neel temperature. The thermodynamic quantities 
calculated are universal functions of the renormalized 
interlayer coupling parameter. The renormalization of 
interlayer coupling parameter turns out to be consider- 
ably temperature dependent. A good agreement with 
experimental data on La2CuC>4 is obtained. The appli- 
cation of the approach used to the case of a ferromagnet 
is discussed. 

75.10.-b, 75.10.Jm, 75.40.Cx 



1. Introduction 

Last time, great interest is paid to properties of quasi- 
two-dimensional (quasi-2D) antiferromagnets in connec- 
tion with the investigations of layered perovskitesu and 
copper-oxide systems, including high-T c superconduc- 
tors. In particular, La2CuC>4 gives one of the best 
known example of a quasi-2D system with small mag- 
netic anisotropy. Unlike 2D systems, quasi-2D ones have 
finite values of magnetic ordering temperature. At small 
interlayer couplings J' the value of magnetic transition 
temperature is small in comparison with the intraplane 
exchange parameter J. There are a number of approx- 
imations which enable to describe the thermodynamics 
of such systems. The standard spin-wave theory (SWT) 
takes into account only the spin-wave excitations which 
exist for quasi-2D systems in a wide temperature range 
up to about J (Refsotj). SWT does not take into ac- 
count the dynamic and kinematic interaction between 
spin waves, which are important at temperatures near 



magnetic phase transition point. By this reason, SWT 
gives too high values of magnetic transition tempera^ 
ture. Recently, the self-consistent spin-wave theoryETcl 
(SSWT) has been proposed which takes into account 
partially the interaction between spin waves. However, 
the value of the Neel temperature in SSWT is still too 
high in comparison with experiment, and the critical be- 
havior is described quite incorrectly. 

To describe the magnetic phase transition we have to 
take into consideration fluctuation (non-spin-wave) cor- 
rections to thermodynamic quantities. It is difficult to 
take into account such corrections in the standard tech- 
nique of the Green's functions because of essentialy non- 
linear character of equations of motion. There exists 
the interpolation approximation by TyablikovO which is 
based on the random-phase decoupling of equations of 
motion for the transverse spin Green's function. This ap- 
proach yields often results which are roughly satisfactory 
from the experimental point of view. At the same time, 
it is difficult to justify and improve such approximations. 

To develop perturbation theory which describes cor- 
rectly the critical behavior, we have to introduce a for- 
mal large parameter in the Heisenberg model. Thus 
the Heisenberg model can be treated as a model with 
a large degeneracy within the 1 /TV-expansion. This ex- 
pansion may be introduced in two different ways. The 
first wayQiia treats the Heisenberg model as a particular 
case (M = 2) of the SU(M) model (i.e. of the model 
with M states per spin degree of freedom at each site) . 
Sinc£ the M — ► oo limit corresponds to SSWT (see, e.g., 
Ref.cl) , at finite M thermodynamics is described in terms 
of the spinrwave picture of excitation spectrum. The 
second wayEH is to consider the Heisenberg model as 
a particular case (N = 3) of the 0(N) model (i.e. of 
the model with TV-component spins) . The limit N — > oo 
gives the quantum spherical model and the large- N case 
corresponds to the fluctuation (non-spin-wave) picture. 
The advantage of the 1/N (or 1/M) expansions over, 
say, the quasiclassical 1/5 expansion is their applicabil- 
ity near the phase transition temperature. 

Since N = 3 and M = 2 are in fact not large, the con- 
vergence of such expansions must be investigated sep- 
arately. For low-dimensional magnets with d = 2 (see 
Refill) and d = 2 + e (Ref.0) the results in the 517 (M) 
model coincide in the zeroth order in 1 /M with those of 
the one- loop RG analysis, and in the first order in 1/M 
with the results of the two- loop RG analysis. In these 
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cases the 1/M-corrections to thermodynamic quantities 
are small. However, quasi-2D systems belong to 3D 
symmetry group so that corresponding 1 /M-corrections 
are not small (see discussion in RefJa) and the series in 
1/M is poorly convergent. Unlike the 1/M-expansion 
in the SU(M) model, the first-order 1/iV-corrections in 
the O(N) model, which were considered in the quan- 
tum 2D caseEJ and in the classical case at an arbitrary 
dimensionality 2 < d < 4 (see, e.g., Ref.Ej), lead to 
results which are close to those obtained by other meth- 
ods. The applicability of the 1/A^-expansion at arbitrary 
dimensionality 2 < d < 4 is important for the investi- 
gation of quasi-2D systems since they demonstrate the 
dimensional crossover from 2D to 3D behavior (see, e.g., 
Ref.El). On the other hand, the renormahzation group 
e-expansion is not appliciablc for d = 2 and d = 3 simul- 
taneously: for e = d — 2 it cannot describe satisfactorily 
the case d = 3 and vice versa, for e — 4 — d the behavior 
at d — > 2 is poor. 

Thus, instead of direct calculation of corrections to 
SSWT, we start in this paper from the quantum spheri- 
cal model, O(oo) and then find the 1/W-corrections. Al- 
though the results in the 0(oo) and SU(oo) models are 
different, it will be shown that already in the first or- 
der in 1/N at low enough temperatures the results in 
the O(N) model are identical to those in SU(oo) (i.e. 
in SSWT). At higher temperatures the results of SSWT 
are modified due to fluctuation corrections. 

The plan of the paper is as follows. In Sect. 2 we re- 
view various approximations in the theory of quasi-2D 
systems, which are based on the spin- wave picture of 
excitation spectrum, and analyze the corresponding ex- 
pressions for the Neel temperature. In Sect. 3 we formu- 
late the O(N) model for the quasi-2D case and the tech- 
nique of the l/A^-expansion, which is a generalization of 
that by Chubukov et all!3 for the 2D case. In Sect. 4 we 
calculate the magnetization, Neel temerature and spin 
correlation function to first order in 1/A^. In Sect. 5 we 
discuss our results and compare them with experimental 
data on La2Cu04. 

2. Spin-wave approximations in the theory of quasi- 
2D antiferromagnets 

We start from the Heisenberg Hamiltonian of a quasi-2D 
antiferromagnet 

H = o ^ JijSiSj (!) 

i'j 

with the exchange interactions Ji^+s = J for S in a plane 
and Ji.i+s = J' for 8 perpendicular to the planes. 



At small values of interlayer coupling J' it is possi- 
ble to derive analytical results for the Neel temperature. 
First we consider the standard spin-wave theory. The 
spectrum of spin waves has the form 

= 5(J 2 - JlY'* (2) 

where J q is the Fourier transforms of the exchange pa- 
rameter 

J q = 2 J(cos q x + cos q y ) + 2 J' cos q z (3) 
The sublattice magnetization is determined by 

q 4 

For small values of J' /J SWT yields different analyti- 
cal expressions for the Neel temperatue in the quantum 
regime (TnccI *C JS) and classical regime (Tn o1 3> JS). 
We have 

T$YJ = 4nJS*x 

f l/ln(T2 ccl /8JJ'S 2 ) ln( J /J') » 2-kS 
\ 1/ ln( Jql/ J') 1 < In (J/ J') < 2nS 

Here qo — tt is a cutoff parameter determined by the 
boundary of the Brillouin zone. Note that for the quan- 
tum case the main contribution to integrals over the 
wavevector comes from the region with q < T , while 
in the classical case the value of TnocI is determined by 
the whole Brillouin zone. 

The spin-wave spectrum in SSWT and the TyablikpM 
approach is renormalized in different ways. SSWTcl cl 
takes into account the interaction between spin waves in 
the simplest self-consistent Born approximation. There 
exist s(peMj generalizations of SSWT on quasi-2D, 
systemstTtJ. We will follow to approach of Refs.liirt£l 
which gives more satisfactory results at small J' / J. The 
spin-wave spectrum in SSWT has the form 

C T = ^7o 2 -tJ /2 (6) 
7 q = 27 (cos q x + cos q v ) + 2-f' cos q z 

where 7 and 7' are the renormalized exchange parame- 
ters which are determined from the self-consistent equa- 
tions 

7/ j - E if cos & coth § + 2 ^( T ) ( ? ) 

y /i7 , = ^ M cos qz coth |S. + 2S(T) (8) 
The sublattice magnetization is given by 
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s = s- 



^^coth^ 



^ 2E n 



2T 



(9) 



V5 = E¥ tanh § 



(15) 



At small values of J' / J we have 



- - T 

S = S 



and has a more complicated form for higher spinsEl. Near 
the Neel temperature TT yields at arbitrary S and any 
ln(T 2 /877'S 2 ) S(JJ') 1 / 2 < T < Jgpace dimensionality d > 2 
4tt 7 S ' ' t ln(79o/7') ^ < T < 



(10) 

where So is the sublattice magnetization in the ground 
state. The quantity 7 varies slowly with temperature in 
the whole region T < TnccI and may be replacecLby its 
zero-temperature value 7(0). According to Refs.ErB, we 
have 



SJ< 



11 



T 



rpTT 

^Ncol 



1/2 



(16) 



where T$ is some function of S, rj/ 2 = 3. Thus, unlike 
SSWT, the critical exponent for the magnetization has 
the standard mean-field value, /3tt = 1/2. For small 
J'/ J, TT yields 



7(0) = 1.1571J, So = 0.3034 



(11) 



rriTT 



for S = 1/2 and 7(0) = J, S = S for S -> 00. The 
second case in ([l0]) may be realized only in the classical 
limit S » 1. One can see from (|l(i| ) that the value of 
the critical exponent for the magnetization is Psw = 1- 
The same critical behavior takes place at an arbitrary 
d > 2. This result is correct only at tL« 2 + e to leading 
order in s (/3 = 1 — 2s, see, e.g., Refs.E!lt!l), and for higher 
dimensionalities (3 <1. 

As follows from (0), the Neel temperature is deter- 
mined by the equation of Ref.L3 (note that some coeffi- 
cients in this paper are incorrect) 



AnJS 2 



(17) 



and 



= 47T7 C SS 



i c = (T Neel /4w lc S 2 )J' 



The result (|17|) is lower than SSWT value 
closed to experimental data (see Sect. 5). On the other 
hand, the result ([H]) coincides with that of the spheri- 
cal mockl.. (which is adequate only in the classical limit 
S — > odlallj) and with the result of the spin-wave ap- 
proximation (||) in the classical regime Tn cc i ^ JS. 
The Tyablikov approximation gives the same result ( |l7| ) 
(with the replacement J — > — J, J' — > — J') for the Curie 
temperature of a ferromagnet (J, J' < 0). This demon- 
strates that near the critical temperature TT does not 
take into account quantum fluctuations which are im- 
(^wrtant for small values of S. Thus we may conclude 
l/ln(T^ eel /87c 7 ^,S 2 ) ln(J/J') > 27rS t nat TT is satisfactory from the practical, but not from 
1/ ln(^ 7c /7^) 1 < ln(J/J') < 27r^icoretical point of view. 

To leading logarithmic accuracy, all the discussed ap- 
proaches give the same value of the Neel temperature. 
However, this accuracy is insufficient to treat experi- 
mental data. In particular, the factor of q$ ~ 10 in 
the classical regime is often not taken into account (see, 
e.g., Ref.tH), although this factor gives an essential con- 
tribution to TNeel- 

To improve the description of the critical region and 
obtain a better appoximation for the Neel temperature 
in the quantum case, it is necessary to take into account 
fluctuation corrections to the spin- wave theory result for 
T N eel © more correctly than in SSWT and TT. To this 
end we use in the next Sections the 1 /TV-expansion in 
the 0(N) model. 



r SSWT 
'Neel 



Here 7 C ~ 7(0) and j' c are the renormalized exchange 
parameters at T = TnccI, The value of 7 £ determined 
from (@) is 



(13) 



in both the quantum and classical regimes. Note that 
the renormalization of the interlayer coupling in ( |l2| ) 
plays a crucial role in lowering the Neel temperature in 
comparison with its SWT value (§) since 7c7 ^./JJ' = 
T Nccl /4irJS 2 « 1. 

In the Tyablikov theory! (TT) the excitation spectrum 
has the form 



E. 



TT 



s(j 2 ~j 2 y 



(14) 



As well as in a ferromagnet, the proportionality of the 
spectrum to S is not quite correct at low temperatures: 
in the antiferromagnet the spin- wave frequency varies as 
T 4 , while the sublattice magnetization as T 2 (see, e.g., 
Ref.O). The equation for S at S = 1/2 reads 



3. The quantum nonlinear sigma-model and O(N) 
model for quasi-2D quantum antiferromagnets 

To describe thermodynamics of quantum antiferromag- 
nets we consider the nonlinear sigma-model which was 
proposed for the one-dimensional Heisenberg model in 
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|r _.In the 2D case this model was applied in 
jO The large value of the correlation length 



Ref.E 
RefsJH 

£ ^> a (a is the lattice parameter in the plane) plays 
the crucial role in the Haldane's mapping of an antifcr- 
romagnetic Heisenberg model (|l|) to quantum nonlinear 
sigma-model (see e.g. RefJ-a). This gives a possibility to 
separate and integrate out the "fast" modes with space 
scale / < A -1 (A satisfies to £ _1 <C A <C cT 1 ) retaining 
"slow" modes with I > A -1 . 

In the quasi-2D case we have £(T < Tn oi) = oo. How- 
ever, at small q we have 



(22) 



J — J q ~ J{aq) 2 + 2J'(1 — cosq z ) 



(18) 



Thus besides the "true" correlation length £, there ex- 
ists also another variable with scaling dimensionality of 
length 



& = l/« 



1/2 



> a 



(19) 



where a = 2 J' / Ja 2 is the interlayer coupling parameter; 
in this paper we consider only the case where a < 1. 
On the scale of order of £j< the regime of fluctuations 
changes from 2D to 3D one. Thus we may use the scale 
£j/ to separate "fast" and "slow" modes in the Haldane's 
mapping. Depending on the value of the imaginary time 
slab thickness 



L T = c/T 



(20) 



(c ~ JSa is the fully renormalized spin- wave velocity), 
three regimes are possible 

(i) L T ~ or, equivalently, T ~ a^^c ~ (JJ') 1/2 . 
This is an analog of the quantum critical regime L T ~ £ 
in the 2D caseLiM 

(ii) a < I T < i.e. a 1 / 2 c <C T -C c which is an 
analog of the renormalized classical regime a -C L r <C £ 
in the 2D case and 

(iii) the classical regime L T «a (i.e. J 5 -C T). 

Since the regime (i) is well described by the standard 
spin- wave theory (or by SSWT), we do not treat the 
thermodynamics at temperatures of order of (J J') 1 / 2 . 
From (|J) and ([l2|) one can see that T Nocl > ( JJ') 1/2 . In 
the regimes (ii) and (iii) implementation of principles of 
finite-size scaling gives 

Tkeel = Ps$(£j> /L T ,£,j> /a) 

were p s ~ JS 2 is the fully renormalized spin stiffness, 
<f>(x,y) is a scaling function with $(oo,oo) = 0. In the 
regime (ii) we have £j> / L T <C O' / a i so that 

T Noel = p s $(l;j,/L T ,ac) = ps^T^/a^c) (21) 

while in the regime (iii) £j/ /£ r ^> £j> /a and 



Note that the results of SWT (§) and SSWT (§J|) for 
the Neel temperature agree with (^Tj) for the quantum 
regime and with (22) for the classical regime. At the 
same time, the result of the Tyablikov approximation 
( |l7| ) satisfies the classical regime scaling form (|2^) for all 
spin values, which confirms the absence of quantum fluc- 
tuations at the critical temperature in this approxima- 
tion. As follows from (^2|), the value of Neel temperature 
in the classical regime depends on fluctuations on a scale 
of order of lattice constant, i.e. is non-universal. There- 
fore in this regime we cannot eliminate "fast" modes by 
Haldane's mapping. Further we will assume that the 
"renormalized classical" regime (ii) takes place. _ . 

We use the same procedure as used by HaldaneE] (see 
full discussion in Ref.Q) to integrate out "fast" modes. 
Thus the partition function has in terms of a path inte- 
gral the form 

il/T 
-y J dr^(a T (T 4 ) 2 
o ' 

l/T ~\ 

~\ S j dT^Jyto-ff^llJ^-l) (23) 
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where <r is a three-component unit-length vector field, i 
is the index of a site, xo is the uniform magnetic sus- 
ceptibility. In the continual limit we reproduce the stan- 
dard three-dimensional quantum nonlinear-sigma model. 
However, in the quasi-2D case the large value of £j< gives 
a possibility to pass to the continual limit only within the 
layers: c j; (r) — > <T^(r,T) where r is a 2D vector, i z is 
the index of a layer. The partition function takes the 
form 



j drj d 2 r^2 



-j(<9 T <T;J 2 



L<-0 



-(y^i.y + t7 <t i - <T, 



1) 



(24) 



where p® — JS 2 is the bare spin stiffness, cq = 
(fls/Xo) is the bare value of the spin- wave velocity. 
Here and hereafter we use the system of units where 
a=l. 

To pass to the O(N) model we replace the three- 
component field (Ti z (r, r) by the A^-component one 
cr™(r,T), m = 1...N. The constraint condition a 2 = 1 
may be taken into account by introducing the slave 
field A,-, (r, r). To calculate the dynamic susceptibility we 
also introduce the external non-uniform time-dependent 
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magnetic field h™(r,r). Then we obtain the partion tunc 
tion of the O(N) model in the form: 



Z[h] = J 



DaDXcxri ( 

^ 2ff 



l/T 



dr / d z r 



J ~E 



spin space, further we will assume W m — aS m N where 
a plays the role of the relative sublattice magnetization 
S/S. Then G NN corresponds to the longitudinal Green's 
^ function, Gi, while other diagonal components (wich are 

— (3,-cri, equal) to the transverse Green's function, G t . At 



I \2 



iK4 1-1) 



-2gh tz (a, 

where g — N/p a s is the coupling constant, ct™ 1 = 
(<7™(r,T)) is the average part of the field c, which is 
supposed to be static and uniform. After integrating 
over a = a — a the partition function takes the form 



Z[h] = 
S eff [X,h} 

2ff 



DXexp(NS e ff[X,h]) 

1 , , a 1 
2^ 

Go (iXa 



2 

(iXa 



In det Go 
-h/p° s ) 



(1 -ff 2 )Sp(iA) 
h/ P ° s ) 



where 



Go = [d 2 T /cl + V 2 



(26) 



(27) 



(28) 



A z cr i2 (r,T) = cr ij8+1 (r,r) - tr ije (r,r) 



Since A enters ( pq) only as a prefactor in the exponent, 
expanding near the saddle point generates a series in 
1/A. At T < TnocI we have the saddle point value iX = 
and a 2 7^ 0. The Green's function of the field cx is defined 
by 



G" m (q,q z 



" ^[0] 

a 2 ^[/i] 



(27T) 2 



dp z 
2ir 



E 



dh m (p,p z ,LUi)dh n (q- p,q z -p z ,cji- n ) 



h=0 



(29) 



; (r,r). 



where h(p,p z ,ui) is the Fourier transform of hi 
Note that only diagonal elements G mm are nonzero, and 
they are proportional to the non-uniform dynamic spin 
susceptibility: 



G mn (q,q z 



£± y r, 



l (q + Q 1 q z +TT,uj)S mn (30) 



where Q = (n, tt) is the wavevector of antiferromagnetic 
structure in the plane; for N = 3 



X 



a/3 



(q,q z ,uj) = ^ 1 



i(qR;+g z i?f) 



(31) 



where Sf are spin operators, a, (3 — x,y,z. Since the 
partition function Z[0] is invariant under rotations in the 



T < TNeei, the value of a is determined by the constraint 
(a 2 ) = 1 which takes the form 



(25) 



a 2 



d 2 k 



dkz 
2tt 



-G r 



l (k, k z ,LJ n ) 
(32) 



We use the relativistic (hard) cutoff w 2 + k 2 < A 2 of 
frequency summations and momentum integrations; in 
this regularization scheme the value of the bare spin wave 
velocity Co is replaced by the fully renormalized one, c, 
which will be putted to be equal to unity except for the 
final results. 

In the limit A — > 00 we may replace in ( |27j ) A by its 
saddle-point value to obtain the "free" Green's function 
(which is the same for transverse and longitudinal com- 
ponents) 



G (k,k z ,u) n ) = [ul + fc 2 + a(l-cosfc z )]" 



(33) 



After evaluation of the integrals and frequency summa- 
tion in ([32J) we obtain the Neel temperature in the limit 

A -> 00 



rpO 



Anp 



N=oo 



iVln(2T 2 ecl /a C 2) 



(34) 



where p^ 00 = N(l/g — l/g c ) is the renormalized spin 
stiffness in zeroth order in 1/N, g c = 2tt 2 /A. To compare 
the result @ with the result of the SSWT we note that 
the value of spin stiffness in SSWT is pf WT = 7 55 
(for 5 = 1/2 this equals to 0.1 76 J which iSrSomewhat 
lower than the result,of two-loop RG analysiscil and nu- 
merical calculations^ 2 ], p s = 0.181J) and the value of the 
spin-wave velocity is c SSWT = V&^/S. Thus we see that 
the value (34) is A times smaller than the corresponding 
SSWT value @ (besides that, in SSWT a is replaced 
by its renormalized value, a^ SWT = 2Y c /'f c < a). r Enr- 
ther we will show that, as well as in the calculation^ of 
the correlation length in the 2D-case in the first order in 
1/N, 

(i) the factor of N in the denominator of (|34|) is to be 
replaced by N — 2 

(ii) p^ =oc and a in ([34]) are to be replaced by their renor- 
malized values, p s and a c 

(iii) terms of order of lnln(2T 2 /a) and unity, which do 
not enter the SSWT result for TnccI, occur in the denom- 
inator of . 

The exact Green's function may be expressed as 
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G mm {k,k Zl uo n ) = [uo 2 n + k 2 + a(l - cos k z ) (35) 



1 = .^E/ ^/lr Go(fc ' fc - £ 



To first order in 1/iV the self-energy £(fc, k z ,uj n ) and the 
function C{k,k z ,oj n ), which describes renoHpalizations where 
owing to the long-range order, are given bytHI 



(27T) 2 

- 3 i?(T,^) + a 2 [1-F(T,^)] 



(43) 



(36) »m J V 



s(fc,fc z ,w„) = — 



A 



<i 2 q /" dq 



2n 



(27T) 



dfc z I(k,k z ,uj m ) 



) 2 J 2tt n(fe, fe z ,c 



(44) 



G (k + q, fc z + g z ,u;„ + u m ) - G (g, q z ,u; m ) 



n(g, g z ,w m ) 



C(k,k z ,u; n ) 

where 

II(g,g z ,w„) 



2a 2 



9 U(k,k z ,ui n ) 



F(T,Xs) 

(37) 



2T 
~N 



E 



and 



ll(g, g z ,w„) + -a 2 G (q,q z ,u n ) (38) 
5 



Ana 2 /gT 



rf 2 k /" rffc z G 2 l (k,k z ,uj m ) 
Wi 27r n(fe,fc z , Wm ) ' 

(45) 



(46) 



n(g, ?2 ,w n )=T^] y y ^-G (p,p z ,uji) 



The calculation of functions /, II for the quasi-2D case 

is presented in Appendx A. 

Thus the functions R and F determine the 1/N- 

corrections to the constraint. The expressions (}43|)-(f45|) 

x Go(p + q,p z + q z , + U) n ). (39) enable one to investigate the magnetization and to cal- 

AT , , ,, ... . /TT^In i • r j. ii m culate the Neel temperature for a quantum quasi-2D an- 

JNote that the quantity C m (p71) has m fact the zeroth . r r n ^ 

_ v| — _ _ r.ifprrnmao'nOT 



order in 1 /N, but the corresponding contribution to the 1 crroma S ne 
constraint is of order of 1/N. The polarization operator 
II(g, q Zl uj n ) determines the longitudinal Green's function 
in the zeroth order in 1/AT 



4. The sublattice magnetization, Neel temperature 
and correlation functions 



<Sf=°°(g,g„0 = 



U(q,q z ,w n ) 



q 2 n(q,q Zl uj n ) + 2a 2 /g 



(40) 



As discussed in the beginning of previous section, wc 
consider the quantum case with a being small enough 
To first order in 1/N the constraint (|32|) takes the form to satisfy the condition ln(2T ] 2 rool /ac 2 ) > 1. The cal- 
culation of the functions R and F at T > a 1 / 2 (i.e. 
(TiP* {JJ') 1 ^ 2 ) is discussed in Appendix B. Neglecting 



l-a s = gT^2 Jj^S j ^Go(k,k z ,u m ) 

~gTj2 / (^ / ^G 2 (M z ,^ m )E(fc, 



the terms of order of 1/ ln(2T^ cel /ac 2 ) we have 

k z ,u m ) mT x T 2T 2 (3 + 2x w )T 4irp s 
R(T, x w ) = — — In ; — In 



2ct 2 T 
N 



E 



rf 2 k I" dk z Gl(k,k z ,u) m ) 

w 2 i ^ m,k z ,u; m ) • 



2ttA 
T 



a AnN 
\n{2T 2 /a) 



NTXn 



Following to@ we introduce the function 

i(k,k Z ,u> m )=Tj2 j y^ G o^^ 



2ttN hi(2T 2 /a)+x w 

8T , 2T 2 , NA 
In In ■ 



w„) (42) 



3tt 3 A 
2T 



■In- 



a 16p 
NA T 



x [G (k + q, /c z + g z , w m + LU n ) 
- G (k, k z ,uj m )] 

and represent the equation (^lj) in the following conve- 
nient form 



and 



3tt 3 A 16p s 4?r 
4irp s 



+ —h(x w ) 



(47) 



In 



AA 



N NTxv n 2 N 16 p 



(48) 
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where the functions ii(x^), J 2 (xj?) are defined in Ap- 
pendix B. After substituting (47), ( |48|) into tha-con- 
straint equation ( |43| ) and using the results of RefO for 
the renormalizcd ground-state sublattice magnetization 
(Jo = a{T = 0) = Sq/S and the spin stiffness p s of a 
quantum 2D antiferromagnet, 



Pa 



A 



1 



In 



AA 



3?r 2 A 16/9, 



=~ [ 1 



32 



In- 



AA 



3tt 2 A 16/9 



(49) 



(50) 



one can see that the sublattice magnetization, being ex- 
pressed in terms of quantum-renormalized p s and Wq, 
still depends on A, i.e. is non-universal. To make the 
sublattice magnetization completely universal we have 
to introduce the renormalized parameter of the inter- 
layer coupling 



1 



In 



NA 



3tt 2 A 16/5 



(51) 



We shall demonstrate below that at low enough tem- 
peratures any regular (nondivergent) terms in the renor- 
malized interlayer coupling parameter are absent, so that 
this is renormalized only due to temperature fluctuations 
at higher T. Being rewritten through the renormalized 
parameters, the constraint equation (f43T) reads 



1 



a 2 



—2 



AT 



A a r A NTx w 



47T/9, 

2 ln(2T 2 /a r ) 
Aln(2T 2 /a r ) + 

1 



(52) 



A n ATa;77 



Note that we have simply replaced a by a r in the terms 
of order of 1 /N in ( |52] ) since this yields an error of order 
ofl/A 2 . 

First we consider the case x-$ 3> 1, or, equivalently, 



NT/Att Ps < a 2 /a 2 



(53) 



Since x-$ is the decreasing function of temperature, this 
inequality is satified at low enough temperatures. In this 
case the integrals I\ {x-s) and I 2 (x&) are of order of 1 /x-^ , 
i.e. are small. Thus to leading (zeroth) order in l/x a the 
constraint equation ( |52| ) coincides with that in the case 
of space dimensionality d = 2 + e (Appendix C) with the 
replacement 1/e — > ln(2/a), which corresponds to the 
limit £ — > with simultaneous cutting the integrals over 
quasimomentum on the scale 1/^j'- Similar to the d = 



2 + e case (Appendix C) we transform the logarithmic 
term in the right-hand side of ( |52"| ) into power and replace 
A -> A - 2. Then we have 



0?/ao) 1/ft 
_ NT 



h(x«)] 
2 . 2T 2 

A )ln 



(54) 



4irp s 
2 ln(2T 2 /a r ) 
~N\n(2T 2 /a r )+x w 



r A 
h{Xa) 



hi 



where, being expressed through the renormalized param- 
eters, 



4:irp s 



{N-2)Ta 2 



(55) 



The "critical exponent" (32, which is the limit of /32+e at 
e — * 0, is given by 



1 



2 A-2 



(56) 



As well as in the d = 2 + e case, two-regimes are 
possible under the condition (|5^) . Consider first the low- 
temperature (spin-wave) region where 



T(N - 2) ln(2T 2 /a)/47rp s < a 2 /a 2 



In this region 



a = (7 



1 - 



T(A-l) 2T 



8irp s 



In- 



a r c 



(57) 



(58) 



At A = 3 we reproduce the result of SSWT (|To|) with 
27Y7 being replaced by a r . The factor A — 1 in ( |58| ) 
has a simple physical meaning: this is the number of 
gapless (Goldstone) modes. We can conclude that in the 
temperature interval (|57]) spin-wave excitations give the 
main contribution to the dependence cf(T). Note that 
the spin-wave result ([58]) can be obtained also from the 
untransformed constraint equation (|52|). 

To demonstrate that in the interval ( |57| ) the experi- 
mentally observed interlayer exchange parameter coin- 
cides with a r we calculate the self-energy E(fc, k z , 0). By 
using (Bra) we get 



X{k,k z ,0) 



8fc 2 



In- 



AA 



3tt 2 A 16p, 
irrespectively of k z . Thus we have 

8 



Gt 1 (k,k z ,0) = k 2 



In 



AA 



= Z~ 



3?r 2 A " 16p, 
[k 2 + a r {l 



(59) 



+ a(l — cosfc z ) 



cos q z 
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We see that the renormalized Green's function differs 
from the bare one by the renormalization factor Z and 
by replacement a — ► a r only. Thus the experimentally 
observed (fully renormalized) interlayer coupling is just 
a r . At higher temperatures the temperature renormal- 
ization of interlayer coupling, which will be calculated 
below, becomes important. 

At intermediate temperatures where 

(N-2)T/4np s < a 2 /al < (N -2)1 'ln(2T 2 / 'a) / '4tt p „ 

(61) 

we have a 2D-like critical behavior of the sublattice mag- 
netization, 



(^o) 1/ft 



T 



9T 2 a 2 

(JV_2)ln + 3m^-2 

a r a 



(62) 



For A = 3 we have (32 = 1, which coincides with the 
critical exponent of SWT and SSWT. However, the term 
with In(c7 2 /<7Q), which is present in (p^), leads to a sig- 
nificant modification of the dependence <r(T) in the tem- 
perature region under consideration in comparison with 
SSWT and leads to a considerable lowering the Neel tem- 
perature. With further approaching the transition point 
the behavior of the order parameter changes to the 3D 
one. 

Consider the temperatures which are very close to 
TnccI, so that a is small enough to satisfy the inequality 
x„ <C 1, i.e. 



a 2 /a 2 « (A - 2)T/4irp s 



(63) 



After expanding ( |52| ) near T = Tn cc i, x-^ = 0, picking out 
the logarithmically divergent parts of I\ (x^) and I2 {xw) 
at small x-^ analytically, and evaluating numerically the 
integrals we have 



1 



T 



Nccl 



a 2 






1 - 


+ - 


8 


7T 


2 N 



1 



In- 



Anp s 



N ~ (A - 2)T Ncc i 
In x w - A 

where Aq = 2.8906/Af. The equation for Tn oc i reads 



(64) 



?Necl = ^Ps 



-3 m- 



(N - 2) In 
4irp s 



2T 2 

z ± Nccl 



- 0.0660 



(65) 



(A - 2)r N cci 

As will be clear below the second term in the denom- 
inator, which is of order of lnln(2T ] 2 recl /a) leads to a 
significant lowering of Neel temperature in comparison 



with SSWT (where only the first term is taken into ac- 
count). We collect separately the logarithmic terms in 
(Q) which comes from the quasimomenta q ^> a 1 / 2 (2D 
regime) and q -C a 1 / 2 (3D regime): 



1 - 



T a 2 

J Neel CT n 



'0 



1 



N (N 



Airp s 



2)7} 



Nccl 



(66) 



1 



8 



TT 2 N 



lnXn 



Unlike the "2D-like" regime, the coefficients at the loga- 
rithms arc different. Transforming the logarithmic terms 
into powers we obtain 



A-Kp s 



(A - 2)T Neel 



03/02-1 



1 



1-Ao 



T 



2;h 



where 



TnccI 

(67) 



(68) 



is the true 3D critical exponent for the magnetization. It 
should be noted that we have not to perform the replace- 
ment A — > N — 2 in ( |68[ ) and other contributions which 
come from essentially three-dimensional integrals. Wc 
get for A = 3 the value fa ~ 0.36. The result @ coip-. 
cides with that of the 1 /N expansion in the (j) 4 modeled 
at d = 3, in agreement with the universality hypothesis. 
The dependence d67j ) is to be compared with that in the 
Tyablikov approximation ([l6|) where (3 = 1/2 and the 
dimensional crossover is absent. 

Consider now the self-energy £(fc, k z ,0) at T = Tkeel- 
At a 1 / 2 <fc< ^*poi tne self-energy has the same form 
as in the 2D caseo with £ being replacing by : 



£(fc,fc z ,0) = fc 2 



r\ In 



AA 
16p s 



1 ln H2T 2 c Ja) 1 
N ln(2fc 2 /a) N 
(69) 



Thus the expression for Green's function reads (G 
G t = Gi) 



G(k,k Z) 0) = ± 



(A - 2)T Neel ^ 2P 



4irp s 



A- 1 
A 



1 — T) In ■ 



a 
AA 



l/(JV-2) 



(70) 



16p s 

a 1 / 2 < fc < T N eel 



At fc <C a 1 / 2 , fc 2 <C 1 the fc-dependence of the Green's 
function changes. After integration and frequency sum- 
mation in ([36]) (which are analogous to the calculation 
of the functions R and F in Appendix B) we have 
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S(fe,fe z ,0) = A x k 2 



A 2 k 2 



(71) 



-g(*" + 5«>m r 



*fc 2 /2' 



Here 



A 2 



77 In 



NA 



16p s 
0.6122/iV 



lining 

N a 



0.4564 

N ' 



(72) 



and 



?? = 8/(3tt 2 A0 



(73) 



is the 3D critical exponent for the asymptotics of the 
correlation function at the phase transition point in the 
first order in 1/N. For A~ = 3 we have ij ~ 0.09. Using 
(H) we find 

G-\k, k„0) = (1 + (fc 2 + y fc 2 ) 1_ " /2 (74) 

fc < a 1/2 ,k z < 1 . 



The quantity 

a c = a(l + A 2 )/(l + Ai) 



(75) 



can be interpreted as the renormalized interlayer cou- 
pling at T = T Nccl . 

Using ( |5l| ) we find the following relation between the 
renormalized coupling parameters at low T and at T = 



1.0686 \ 
N J 



(A" - 2)T r 



Nccl 



Airp s 



l/(N-2) 



(76) 



When deriving (|76j) we have transformed the term with 
lnln(2T 2 /a) into a power and then replaced A" — * N — 2 
in the exponent. As well as in SSWT (see Sect. 2), the 
renormalized interlayer coupling at Tn oc i is lower than 
the low-temperature one, but the concrete expression at 
N = 3 is shghtly different from these in SSWT. 

Using (|7^) we get the following equation for TnccI in 
terms of a r 



'Nccl 



Airp, 
+2 In 



2T 2 

(N - 2) In 



(77) 



Anp s 



(N - 2)T Nccl 



1.0117 



where c is the fully-renormalized spin-wave velocity; in 
SSWT we have c = V8j{0)S (see Sect.2). For N = 3 we 
have 



Nccl 



47T/9 S / In 



5.5005 



(47T/9 S 



(78) 



which is similar to the result of the Tyablikov approxi- 
mation ([It]), but the bare value of a is replaced by its 
renormalized value at the critical temperature ([76]) and 
p s is also replaced by its renormalized value. Besides 
that, the result ([78]) does not violate the scaling form 

Finally, we consider the spin correlation function 

w/t^*' (79) 

— 77 

x W ,kImr(M^)^ 

m 

at T = TnccI- For N = 3 we have 

S(R,R Z ) = I < S,(r)S l+i? , 2 (r + R) > | 

The static approximation is sufficient to determine the 
asymptotics of the_cprrelation function. Using (|30|), ( p9| ) 
we derive (cf. Ref.B) 



2 

S(R,R Z )- TS ° 



Ps 



1 + r] In 

7T 

d 2 k 



NA 



{2-Kf 



16p 

2tt 



(80) 



G(fc,fc z ,0)e 



One can see that at R 2 a+R 2 3> 1 the asymptotics of the 
integral in ( ]so|) is determined by the region k <C a 1 / 2 and 
k z -C 1 where G(fc, fc z , 0) is calculated above (see (74)). 
Substituting ( f74| ) into (p0|) we have 



1 TNecl'S'o 
47T p s 



(1 - A, - V C) 



1/2 



1 + rj In 



A^A 



16p s 



(81) 



where 71 = (R 2 +2R 2 z /a c ) 1 / 2 and C ~ 0.5772 is the Euler 
constant. Using the value of A\ ( (72]) and transforming 
the term with 1 /N In ln(2T 2 /o;) into power we obtain the 
final result for the spin correlation function at T = TN ee i 
and TlaJ 2 > 1 : 



S(TZ) = 



TnccI'S'o 

4irp s 



{N - 2)T Ncol 



-, l/(N-2) 



Anps 



xil-Ax-^C) 



1/2 



(82) 
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where A x = 0.4564/iV. Thus S(K) enables one to de- 
termine the value of a c . As one should expect, being 
rewritten through renormalized parameters p s and a c , 
S(1Z) does not contain the cutoff parameter A and is 
thereby completely universal. 

At 1< R < a J 2 we derive from @, @ the leading 
term of asymptotics of the correlation function within a 
plane (N = 3) 



S(R,0) = f- 



T 



hr 



~R 2 



(83) 



Thus we have in this case a logarithmic decrease of the 
correlation function, as well as in the 2D case at 1 <C 



5. Discussion and conclusions 

In the above treatment we analyzed the sublattice mag- 
netization S of a quasi-2D quantum antiferromagnet 
(T N eei < JS). At temperatures T < (JJ') 1/2 the behav- 
ior S(T) is satisfactory decribed by the standard spin- 
wave theory For T ^> (JJ 7 ) 1 / 2 we have obtained the 
equation (|52j) which determines S to first order in the 
formal small parameter X/N. We have three temperature 
intervals (the boundaries of the intervals are presented 
for N = 3) 

(i) the case of low temperatures 



T\u(2T 2 /JJ')/{4ttJS 2 ) < S /S 



[So = S(T = 0); So = S - 0.196 for the square lattice] 
where the results of SSWT are reproduced, 
(ii) the case of the intermediate temperatures (|6l|), or 
equivalently 



S 2 /S 2 < TM2T 2 /JJ')/(AttJS 2 ) 
1 - T/T Ncc i » (1 - A a )(T/47rJS 2 )^ 2 



(1 — Aq ~ 0.0365), where a 2D-like critical behavior, 
which is similar to that in SSWT, takes place. However, 
the corrections to SSWT modify considerably the numer- 
ical factors, so that the Neel temperature is considerably 
lowered 

(iii) the vicinity of the Neel temperature (|63|), or 



1 - T/T Nccl « (1 - A )(T/4ttJS 2 )^ 2 



obtained within the first order in 1 /N, since the equation 
( p2|) is transformed in different ways in these regions to 
derive the results (|54|) and (Q) respectively. Note, that 
in the region (ii) the "2D-like" behavior of the system 
enables one to calculate corrections to SSWT in a regu- 
lar way, e.g., by using the l/7V-expansion in the CP 1 ^^ 1 
model. 

We have also derived the expressions for the mag- 
netic transition temperature ( |65|) , ( |77| ) which contain the 
renormalized quantities a r . c = 2^' r c /j where c are the 
experimentally observable (renormalized) interlayer ex- 
change parameters at low temperatures and T = T^ ee i 
respectively, and 7 ~ 1.1571J is the value of renormal- 
ized intralayer coupling parameter which is weakly tem- 
perature dependent. Therefore these expressions have 
an universal form, in agreement with the scaling anal- 
ysis. Unlike the corresponding results of the spin-wave 
approaches (see Sect. 2), they contain the terms of order 
of In ln(T^ ocl / J J') and unity, which are formally small as 
compared to the leading term of order of ln(T^ eel / J J'). 
However, the lnln-terms result in a significant lowering 
of TNeei in comparison with the SSWT value ( |Io|) at not 
too large ln(J/J'). The regular terms yield small correc- 
tions only, so that one may expect that the higher-order 
terms in 1/ ln(T^ ecl / J J') may be neglected, j— . 

The experimental temperature dependences of the 
sublattice magnetization in La2Cu04 is shown in Fig.l. 
For comparison, the results of spin- wave approximations 
(SWT, SSWT and the Tyablikov theory, see Sect. 2) and 
the result of 1/N expansion are also presented. The 
renormalized value of the in-plane exchange parameter 
7 ~ 1850-ftf can be found from the experimental dataEil 
and the value 7^/7 = 5 • 10~ 4 was chosen from the best 
fit to SWT at low temperatures T < 100 K). The experi- 
mental results for 7V7 are not reliable, and it is difficult 
to compare our value of-^y'/7 with experiment. For ex- 
ample, the result of Ref.tj, 7V7 — 5- 10 -5 , is by an order 
lower than that found from the fit in the spin-wave re- 
gion. It is also important that the value of 7V7 has an 
appreciable temperature dependence because of renor- 
malizations. In particular, we have from ( |76| ) for above 
parameters a c /a r — Y c /y r — 0.13. Thus experiments at 
different temperatures may give different results. 

One can see that SWT and SSWT yield satisfac- 
tory results for T < 0.6TNeei and T < 0.8Tn oc i respec- 
tively. At higher temperatures the sublattice magne- 
tization in SWT and SSWT is still linear in temper- 
ature, so that the critical exponent is /?sw = 1, in- 

(86) 

stead of the experimental one, /3 oxp ~ 0.33. Besides that, 



(84) 



(85) 



where we obtain the critical behaviour S ~ (TnocI— T)^ 3 , 
3 ~ 0.36. 

The detailed description of the temperature region be- 

2 2 

tween (ii) and (iii), where S /S ~ T/4:TTp S: cannot be 



cxp 

both theories give large values of the Neel temperature 
TjSWT = 672K,TSSWT = 537K Thig fact ig often nQt 

taken into account when treating experimental data. At 
the same time, TT gives the value T^? el = 454 K which 
is much lower than those in SWT and SSWT and the 
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magnetization critical exponent /3tt = 1/2. Thus the 
Tyablikov approximation seems to describe the experi- 
mental data more satisfactorily. However, this approxi- 
mation may be justified in fact only in the case of "clas- 
sical" magnets with T Nco i 3> JS. Besides that, TT has 
a number of drawbacks: mean-field values of critical ex- 
ponents, absence of crossover from 2D-like to 3D behav- 
ior of magnetization, neglect of quantum effects at high 
temperatures (in particular, TnccI = ?Curie for the same 
\J\)- 

The result of the 1/iV-expansion to first order in 1/N 
( |65j ) is TnccI = 345 K which is considerably lower than 
in the Tyablikov approximation and is in a good agree- 
ment with the experimental value, T^ooi = 325K. The 
spin-wave region extends up to 300A", while the 2D-like 
region from 320K to about 340K ; the critical 3D region 
is narrow (about IK) . The results of the numerical solu- 
tion of equation ( |54[ ) in the temperature regions (i) and 
(ii) and the dependence ([34]) in the region (iii) turn out 
to be smoothly joined. One can see also that the result 
of the 1/iV-expansion is most close to the experimental 
data and demonstrates a correct critical behavior. One 
may assume that higher-order 1/iV-corrections will give 
a precise description of the experimental situation. Thus 
we may conclude that using the 1/iV-expansion in the 
O(N) model improves considerably the results of stan- 
dard spin- wave approximations in the Heisenberg model. 

Recent experiments demonstrate existence of a gap 
for the out-of-plane spin- wave excitationso in La2Cu04, 
which is assumed to be determined by the easy-plane 
anisotropy. The possible role of-, easy-axis anisotropy 
was also discussed, see, e.g., Ref.EZl. Therefore an exten- 
sion of the present approach to 2D systems with a weak 
anisotropy is of interest. The results may be expected 
to be similar to those in the quasi-2D case, since SSWT 
gives similar descriptions of both_the types of magnets 
with small ordering temperatures, 

The case of "classical" spins cannot be treated consis- 
tently in the continual limit since in this case the natural 
upper limit cutoff parameter (which is the temperature 
in the quantum case) is absent, and the integrals are 
determined by the whole Brillouin zone. Therefore the 
continual models may be used to calculate the critical 
exponents, but not the temperature dependence of mag- 
netization in a broad interval and the Neel temperature. 

It would be also interesting to perform similar calcu- 
lations of thermodynamic properties for a ferromagnet. 
The results should coincide with those for a antifcrro- 
magnet only in the classical case. Unfortunately, the 
nonlinear-sigma model for ferromagnet has the Berry 
phase term A{cr)dcr / dr in the action (A is the vector 
potential of unit magnetic monopole), see, e.g., Ref.cl. 
This term cannot be eliminated in the quantum case 
and prevents constructing the l/TV-expansion. For sin- 



gular contributions, the results for a quantum ferromag- 
net (Tcuric <C JS) differ from those for an antiferromag- 
net by the replacement 



ln(T 2 /85 2 77 ; J -» \n{T/j' rc S) 



(87) 



(as well as in SSWT, see Rcf.0). Taking into account 
only such terms, the expression for the Curie tempera- 
ture has the form 



Curie 



= 47TjO s 



(N - 2) In 



+ 31n 



27rp s 



(N - 2)T Curie 



+ 0(1) 



or, in terms of the renormalized exchange parameter at 
the Curie temperature, 



T. 



Curie 



47rp s 



(N - 2) In 



I'cS 



2 In 



2irp s 



(N - 2)T< 



where 



l 7 7 r 



(n - 2)r Curi , 



l/(N-2) 



C uric 

(89) 



(90) 



and Ay ~ 1. One may expect that, as well as in (|65|), ( |77j ) 
the non-singular terms will influence weakly the value of 
the ordering point. These regular contributions may be 
calculated for a ferromasnet within the l/7V-expansion 
in the SU(N) model (effl). However, as discussed in the 
Introduction, this expansion gives poor results at not too 
large N for d not too close to 2, so that only the "2D-likc" 
region can be described satisfactory. The description 
of the 3D critical behavior of a quasi-2D ferromagnet 
requires other methods. 

We are grateful to M.I.Katsnelson for stimulating dis- 
cussions. 



Appendix A. Analytical results for the functions 

U(q,q z ,w n ) and I(q,q x ,u„) 

Here we present a list of results for the polarization oper- 
ator H(q, q z ,u> n ) ( pS9| ) and the function I(q, q z ,uj n ) deter- 
mined by ( p2| ) at <?/ 2 <C T , and the asymptotic forms 
of these functions. 

Due to inequality a 1 / 2 <C T, the possible values of 
q, u) n may be divided into two regions. The first region 
is uj n = 0, q <C T, while in the second region q 2 +u; 2 ^> a, 
i.e. either uj n — and q 3> a 1 / 2 or uj n ^ at arbitrary 
<?• 

It may be checked that at q <C T and oj n = the main 
contribution to n comes from the term with uj m = 0. 



Oil) 
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After integrating over hJt z with the use of the Feynman 



identity (see, e.g., Refjfij) we get 
1 

dz 



2irq 4 



In 



2q 2 3 + cos q z 



, a < q 2 < T 2 
(98) 



y/z{l-z) y/z(l- z)q i + 2a(q,q z ) 

q < Tl(q,q z ,u 



In 1(]9(l Region with g 2 + u 2 ^ 3> a we have 



Tq 2 



% 2 +^ ) 2 



■In 



where a(g, g z ) = afg 2 + a(l — cos g 2 )]. At laxge q 3> a 1 / 2 
the function H(q,q z ,0) has a 2D form (cf.E3) 



2^ + ^)3 

, T u)% — 3q 2 — (q 2 + u> 2 ) cos q z 



(99) 



4~ 



(<Z 2 +w 2 ) 3 



T 2a 2 

n(g,^,0)~ jj—jln-^-, a«q 2 «T. 



27r g 2 



a 



(92) 



+I qu (g,w n ), q 2 +cj 2 > 



with 



In the opposite limit the form of function H(q,q z ,0) 
changes from 2D to 3D one: 



T 



Iqu(q,Un) = — 77^ 

* ^ [( W 2 + ^ + 2WmW „)2 + ^j8/» 



(100) 



n(9,fe,o) 



V2a(g, g z 



<f < a. 



(93) 



x Arctanh • 



<? 2 + 2w m u;„ 



Consider now the case q 2 + ui 2 3> a. Picking out the 
terms with m = and m — —n (if n 7^ 0) we have 



+ q 2 + 2u m uj n ) 2 + \q 2 u 2 m 
q 2 + 2m m uj n + LP 2 , 



--E 

Air ^—^ 



IL(q,q z ,u„ 



_T m 2(9_+ay_ + , Wn ),In(9*p ultraviolet limit g 2 + cu 2 » T 2 the following 

27r(q- + ujfj ' 



q 2 a 



q 2 + LU 2 t > a 



asymptotics takes place: 
Tq 2 



where the quantum contribution Il qu is given by 



■In- 



n qu (g,^n) = — 2J 



1 



m/0 

x Arc tanh 



^/(u; 2 + <? 2 + 2w m w„) 2 + 4g 2 w 2 
u> 2 + q 2 + 2uj m iv 



Iqu(q,U>n) 
(95) ' 7T ( g 2 +CJ 2)3' 



qT 



n{q 2 + ujI) 3 q 2 + uj 2 
T u 2 n -q 2 



(101) 

7 2 + UJ 2 > T 2 



V ( w n + <7 2 + 2w m w„) 2 + 4g 2 w 2 Appendix B Calculation of 1/iV-corrections to the 



In all the further calculations we will need only the 
asymptotic form of II qu (q, u n ) for q 2 + w 2 ^> T 2 . In this 
limit we find 



n qu (<?,w n ) 



In- 



ir{q 2 +u 2 ) q 2 +uj 2 

+ ; 1 g 2 +^ 2 »T 2 
Sv^T^ 2 ' 



(96) 



For oj„ = and g<T we obtain by analogy with the 
calculation of H(q, q z ,0) the result 



i{q, q z ,o) 



T 



1 



47r q 2 + a(l — cosq z 
1 



constraint at T < TNeei- 

Consider briefly the calculation of the functions R ( [44| ) 
and f ( [45| ) which determine, according to (43), the cor- 
rections to the constraint to first order in 1/N. First 
we introduce intermediate cutoff parameters C and C 
determined by a 1/2 < C < T < 27rC" < A and divide 
the region of summation and integration q 2 + ui 2 < A 2 
into four regions: 

1. u„ = 0, q < C. 

2. LU n = 0, C < q < 2nC 

3. uo n f 0, q 2 + OJ 2 < 2nC 

4. 2ttC" < tfUr uj 2 < A 

Further we denote the contributions from i-th region to 
R and F as Ri(T, x w ) and F,(T, x w ) 



r dz q z(l — z) + az{\ — cos q z ) i n the/fipt region we can use the expressions for the 
J y/z(l-z) [q*z(l -z) + 2a(q,q z )f 2 functions tl(q,q z ,0) andJ(g, q^, 0) at q 2 « T @, © 



and their asymptotics (B2j), (B3|), (Bq). Then 



and its asymptotic form 
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Rl (T, x„) 



T , 2C 2 
In 

2ttN a 



x In 



(3 + 2x w )T 
AttN 

\n(2C 2 /a)+x w T 



(102) 



Contributions from other three regions are also easily 
calculated: 



4tvN 



h(x w ) 



1 



4irp s /NT 



where 



F 2 (T,^)^-ln ln(2C2/a)+ ^ 
F 3 (T, X7 ?) = 0(l/ln(2T* e Ja)) 



h (x w ) = -f- / qdq I 



A' 



2tt 



1 



1 — cos q z 



8 , NA 
m ■ 



l&p s 



(110) 
(111) 
(112) 



I(q,q z ,0)-U(q,q z ,0) 3 + 2x w 6(q 2 - 1/fihe 



Summing up all (i = 1...4) we find the result (E8|) of 



n(«,g«,o) 



2g 2 ln(2g 2 ) 



mam text. 



#(:r) is the step function. In the second and third regions 
we use the expressions for the functions U(q,q z ,0) and 
l(q,q z ,0) at q 2 > a, <M) and ©: 



2 ^m 



In 



, T 3 



2ttC" ln(2T 2 /a) 



ln(2C 2 /a)+x 5 



T ln(2T 2 /a) 



R 3 {T,x w ) 



T \n(2T 2 /a) 



ttN \n{2T 2 /a)+Xe 
where we have used the identity 
ln(2T 2 /a) + x w = Aitp s /NT 



4C" 
3T 



1 , 2ttC 
2- ln — 



(106) 



Appendix C. The order parameter and transition 
temperature at d = 2 + e 

In this Appendix we consider the calculation of the sub- 
lattice magnetization to first order in 1/N in the space 
with the dimensionality d = 2 + e. We will be interested 
in the terms of the leading order in e at not too small 
wMi.'pferatures T ^> Je -1 / £ (which is an analog of the 
renormalized classical regime in the 2D case), so that 
c^mjjgjthe contributions with zero Matsubara frequences 
will be taken into account. Consider first the results for 
the functions II and /. Evaluating the integrals in ( |39| ) 
and (^) at an-arbitrary space dimensionality 2 < d < 4 
(see, e.g., Ref.Ea for the procedure of calculation of such 
integrals) we have 



which is satisfied in the zeroth order in 1/N. In the 
fourth region we obtain 



Ri(T,Xa) 



8T , 2T 2 , NA 
: In In ■ 



3lT 3 N 



2T 



In 



a I6p s 
NA 4C" 



n(g,o) 



TK d A d 
q i-d 

TK d A d (3-d) 

n 6-d 



hi(2T 2 /a) 



3ir 3 N 16p s 37rA^ln(2T 2 /a) 



where we hav e us ed the asymptotic forms (g6J), (101) and 
the identity ( |l06| ) . After collecting all R z (i = 1...4) the 
intermediate cutoff parameters C, C are canceled and 
we find the result (|47|) of the main text. 

Analogously, we obtain the contribution from the first 
region to F in the form 



(WK^re q is the d-dimensional vector, 

r(d/2)r(2-rf/2)r 2 (d/2- 



2T(d-2) 



(113) 
(114) 

(115) 



2 d -\ d/2 T(d/2) 



T(x) is the Euler gamma-function. At d = 2 + e we find 
to leading order in e 



Fi(T,Xff) 
where 



1 ]n(2C 2 /a) 
N 11 xV 



OO 7T 



+ h{x w ) (108) 



1 



' me) 

eq z e 2tt 



_(q 2 + l 
1 9(q 2 



cos q : 
1/2) 



n(q,g z ,0) 2q 2 ln(2q 2 ) + x w 



(109) 



U(q,0) =q 2 I(q,0) 



The constraint equation to first order in 1/N ( |13| ) takes 
the form 



gT l +*K 2 2 
e [ N 

T^/e + a 2 /g 



i-^m 

N 



(117) 
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Using the identity 

gT 1+s / 



a 2 



1 



(118) 



which is satisfied in the zeroth or der in 1/N and trans- 
forming the logarithmic term in (117) into a power, we 
obtain 



[3 2+e = (l + l/N)/2 + (D(l/N 2 ,e) 
The Neel temperature is determined by 



rpl+e 
J Ncol 



6 2irp s e 
gK 2 (l - 2/N) = N-2 



(119) 
(120) 



(121) 



This result coincides with the result of the RG analysis^ 
The RG result for the critical exponent /3 reads 



1 



N 



0(s) 



(122) 



Thus one have to replace N — » N — 2 in (12C). Such 



a replacement is analogous to this in the renormalized 
classical regime of Ref.E3 and may be justified by the 
calculations of terms of order of 1 /N 2 —which we did not 
carry out. As demonstrated in Ref.Ej by calculations 
of analogous contributions up to 1/N 2 , this replacement 
should be indeed performed. Since the denominator in 
(121) is of order of N, the replacement N — > N — 2 occurs 
already in the first-order expression for the transition 



temperature (121^ 

According to (119), two regimes are possible in the 
temperature dependence of the order parameter. At 
Je~ x l e <C T <C TnccI we have the spin- wave behavior 



a = l- 



(N - l)T 1+e 
4irp s e 



(123) 



For N = 3 this result is analogous to the quasi-2D case 
result (|l0|) in the quantum spin case. At 1 — T/TnocI -C 1 
the temperature dependence of the sublattice magneti- 
zation changes from the linear one to the power behavior 
with the critical exponent P 2 +e- 

Two above temperature regimes correspond to differ- 
ent pictures of the excitation spectrum. In the low- 
temperature regime T Tkcoi we have from (|4C| ) at 
quasimomenta q < T (only such q give a contribution to 
thermodynamic quantities) the zeroth order longitudinal 
Green's function 



G?=°°(q,0) = »-Il(q,0) = 



gTK 2 

eq 2 ~ £ 



(124) 



which corresponds to spin-wave excitations. Near the 
phase transition point we have at an arbitrary q (except 



for the exponentially-narrow hydrodynamic region q < 
(2a 2 /<?) 1/£ ) 



G 



N=oo 



(9,0) 



(125) 



which corresponds to critical (non-spin-wave) excita- 
tions. 
Figure caption 

Fig.l. The theoretical temperature dependences of 
the relative sublattice magnetization S/Sq from different 
spin-wave approximations and from the l/TV-expansion 
in the O(N) model (equations (54) and.(p4)), and the 
experimental points for La2Cu04(Rcf.Ej). 
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